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Direct photon production in minimum bias d+Cu and d+Au and central Cu+Cu and 
Au+Au collisions at center of mass energies y/s = 62.4 GeV and 200GeV at RHIC is 
systematically investigated. We study the jet quenching effect, the medium-induced pho- 
ton bremsstrahlung and jet-photon conversion in the hot QGP. We account for known 
cold nuclear matter effects, such as the isospin effect, the Cronin effect, shadowing and 
cold nuclear matter energy loss. It is shown that at high pT the nuclear modification 
factor for direct photons R\j^^{pt) < 1 is dominated by cold nuclear matter effects and 
there is no evidence for large cross section amplification due to medium-induced photon 
bremsstrahlung and jet-photon conversion in the medium. Comparison of numerical sim- 
ulations to experimental data also rules out large Cronin enhancement and incoherent 
photon emission in the QGP but the error bars in the current experimental data cannot 
provide further constraints on the magnitudes of other nuclear matter effects. 
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1. Introduction 

In relativistic heavy-ion collisions a large amount of energy is deposited in the 
interaction region over a very short time interval t ~ 2Ra/"1- An excited state of 
matter, the quark-gluon plasma (QGP), is expected to be formed and its properties 
are subject to intense theoretical and experimental in vesti gation. From SPS to RHIC 
and to the LHC, with increasing ^/s, hard probes have become increasingly 
more important as tomographic tools in the study of the QGP. Their perturbative 
QCD (pQCD) calculation is based on the factorization theorem, which separates the 
hard partonic part from the soft, non-perturbative part. For a physics observable 
dF/dy(Pprp in p+p collisions we have: 

dF 

« <Pa/Ai^A, ^J'f)(|)b/Bi^B, fJ-f) 



dyd'^ pi 

^ dyd^pr (1) 
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Here, dF"'^'''^'^'^^ /dycPpT represents the partonic cross section, while 4>a/A{^A, fJ-f) 
and Dii/c{z, fJ^f) are the non -perturbative parton distribution functions (PDF) and 
parton fragmentation functions (FF), respectively. In nucleus- nucleus (A+A) col- 
lisions, the result of final-state interactions can often be represented as effective 
modifications of the FFs. Similarly, cold nuclear matter (CNM) effects resulting 
from the initial-state interations in cold nuclei can be absorbed in effective mod- 
ifications of the PDFs. The latter will also be present in A-l-A collisions and a 
robust calculation of QGP signatures needs reliable understanding of initial-state 
interactions. 

Since the photon couples to the partons in the collision region only through elec- 
tromagnetic interactions (aem ^ Q^s) its mean free path is very large. Hard photons 
leave the medium without rescattering and their production can provide valuable 
information to help constrain initial-state CNM effects. The study of 7 production 
in A+A has so far been centered on its QGP-specific components, in particular 
medium-induced photon emission in the QGP and jet-photon conversion The 
predicted consequential large enhancement of the direct photon production cross 
section in central heavy-ion collisions appears to be incompatible with recent data^l. 
Our paper ^ addresses the need for a systematic study of direct photon production 
in both p-|-A and A-|-A reactions by taking consistently into account all relevant 
hot and cold nuclear matter effects. We concentrate on the intermediate and large 
transverse momentum regions {pT > 2 GeV) in minimum bias d-f Cu and d+Au 
and central Cu-t-Cu and Au-|-Au heavy ion collisions at RHIC with center of mass 
energies of y/s = 62.4 GeV and 200 GeV. We also provide the theoretical derivation 
of the medium-induced radiative 7 spectrum for a quark emerging out of a large 
process. 

2. Direct photon production in p+p collisions 

In "elementary" p-l-p collisions direct photons can be produced, at leading-order in 
perturbation theory, by quark-antiquark annihilation {q + q 73) and Compton 
scattering (q+g — > jq), see LHS of Fig.[TJ One should also include the contribution of 
the bremsstrahlung/fragmentation process illustrated in the RHS of Fig.[TJ Naively, 
it seems that 7 radiation gives a higher-order contribution {aemOtl) when compared 
to annihilation and Compton scattering (aemds)- However, noticing the logarithmic 
growth with the hard scale of the photon fragmentation function D^^c{z,Q'^), 
parametrically we have: 

a^iQ) cx ln"^(^), D^^/ciz, Q^) cx In(^) 
aemal{Q^)D^/^{z, Q^) cx ae™as(Q^). (2) 

''The convenience of a particular mathematical representation of nuclear matter effects on the 
observable cross section should not be confused with genuine universal modification (which may 
also be present) of the parton distribution and fragmentation functions. 
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Fig. 1. Direct photon production from the q + q annihilation process + Compton scattering (left 
panel) and bremsstrahlung/fragmentation (right panel). 




Pt [GeV] 



Fig. 2. Cross sections for direct 7 production in P+P collisions at y^i = 62.4 GeV and 200 GeV. 
Data at = 200 GeV is taken from PHENIXl^. Insert gives the ratio of bremsstrahlung photons 
to all direct photons. 



Hence the contribution of the bremsstrahlung process is effectively the same order 
as the prompt production processes and should be taken into account even in LO 
calculations. 

In our leading-order pQCD model, the cross section of direct photon production 
in p+p collision is given by Eq. ([T]) , where we also allow for possible deviations from 
collincarity for the incoming partons: 

(X j SliadPU f{ka)f{kb){---) . (3) 

Further technical details are given in ^ and we emphasize that any overall scale 
factor (e.g. the process-dependent NLO K-factor) will cancel when we calculate the 
nuclear modification ratio Rab{pt)- 

In Fig. [2] we show our numerical results for the differential cross sections for direct 
photons in p-l-p collisions at ^/s = 62.4 GeV and ^/s = 200 GeV. The power-law 
Pt dependence of the data at 200 GeV, measured by PHENIX is well described 
by the pQCD calculation. The insert in Fig. [2] gives the fraction of bremsstrahlung 
(fragmentation) photons to all direct photons. We note that at very high pT the 
fragmentation 7 yield 25% — 30% of all direct photons in p-|-p collisions. 
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3. Direct photon production in heavy-ion collisions: hot nuclear 
matter effects 

We first discuss final-state QGP-induced modifications to da"^ /dycpTpT- 

• Jet quenching in the QGP 

It is well established that when an energetic parton propagates in a 
hot /dense nuclear medium, it will suffer multiple scattering and lose a frac- 
tion of its energy via induced gluon radiation This jet quenching effect 
will reduce the contribution of the fragmentation photons due to the at- 
tenuation of the flux of energetic quarks. We calculate the parton energy 
loss using the GLV formalism and represent the effective fragmentation 
functions into photons as follows: 

D^,,{z) ^ \t P{e) -^D^/c (y^) . (4) 

Here, P(e) is the probability distribution of the fractional jet energy loss 
e = l:^E/E^. 

• Medium-induced photon radiation 

In the hot QGP, a propagating jet may also radiate photons due to mul- 
tiple scattering in the medium Because there is no gauge-boson self- 
interactions in 7 bremsstrahlung, one may expect that the dl^ /duj evalua- 
tion is a trivial simplification of the gluon intensity calculation. This naive 
expectation is not true. Consider the radiative amplitude for single scatter- 
ing of a fast on-shell quark (there are 3 diagrams in QCD, while 2 diagrams 
in QED as shown in Fig. [S]): 

M'Uk) cx 2^g,e^ ■ (^^ - j^-Sji^ e^^^" ^ . (5) 

In the QED hmit T'^ ^ 1 and we obtain Ml^j^{k) -> 0. Thus, a theo- 
retical approach optimized for medium-induced gluon radiation cannot be 
directly applied to photon bremsstrahlung in the medium and one should 
be careful i n id entifying the correct kinematic approximations in evaluating 
7 emission ISEl 




Fig. 3. Single-Born diagrams for medium-induced 7 emission. 
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Meaningful direct photon phenomenology requires theoretical advances in 
understanding the medium-induced 7 radiation off of quarks produced in- 
side the plasma in large processes. We derive this contribution to the 
7 spectrum using the Reaction Operator approach EE. For the single-Born 
scattering diagrams, illustrated in Fig. [3l we obtain: El 

M.„,„,(,), = e(i-a_i_£i),..f*-, (e) 

where the collisional part of the amplitude is not shown. The contribution 
of Double-Born diagram is found to be negligible: M^^^{k) « . Extend- 
ing the above calculations to higher orders in the correlation between the 
multiple scattering centers we demonstrate that contributions to dP /(Lj 
vanish beyond second order in opacity El 

dN'^{k) _aemf f dAzi /•^2q 1 '^^cr"' 



X 



dk+d^k^ TT^\JXg{zi)J ^^a^^d-^q^i 

]Mradi{m' + 2M:,,i{l})Mradm) COs(fc- Az+)] 



-|- corrections . (7) 

Here Az^ = zf — z^_i, and Xq is the quark mean free path. We define the 
inverse photon formation time as ry^ = k^/(2a') « ^2k~ . There are two 
well-defined limits in Eq. ([7]): when T^^Xq 3> 1, the term cx coa{k^ Az^) 
vanishes due to oscillation and we recover incoherent 7 emission. In the 
opposite case, ry^Xq < 1, interference between the vacuum and medium- 
induced photons becomes critical. 




Au+Au at b-3 fm, s - 200 GeV 



yfrom u, d quark - incoherent Btt 
ytrom c quark ■ incoherent Btt 
ytrom b quark - incoherent Btt 
yfrom u, d quarks - coherent FS(1) 
yfrom c quarks - coherent FS(t) 
yfrom b quarks ■ coherent FS(1) 



0.01 
X = k* / E* 
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Fig. 4. Scaled medium-induced photon number spectrum versus x = / E'^ for Eq = 100 GeV 
light, charm and bottom quarks in central Au+Au collisions at y/s = 200 GeV. Solid curves 
represent the incoherent Bethe-Heitler (BH) photon radiation in the medium, and dashed curves 
stand for the coherent final-state (FS) photon bremsstrahlung. 



In Fig. m we show numerical results for the scaled medium-induced photon 
number spectrum dN'^ /dx = {e/eq)'^dN'^ /dx, x — k^ /E'^ for a quark jet 
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propagating outwards from the center of the medium created in Au+Au col- 
Usions at RHIC with the impact parameter 6 = 3 fm . It is apparent that in- 
terference effect strongly suppresses stimulated photon emission when com- 
pared to the incoherent scenario. This reduction in emission strength, which 
is relevant for photon production in heavy ion reactions, differs significantly 
from the argued small LPM ef fect for on-shell quarks, see e.g. 121101 jj^ q^j. 
detailed numerical simulations '^1^ we also find that the average number of 
jet interactions in the medium is never large, (n) = L/Xq = 2 — 3. 

• Jet-photon conversion 

In the hot /den se medium, a possible new source of direct 7 is jet-photon 
conversion I^ISI^ gee e.g. Fig. [5l Making the approximation that p-y ~ Pc 
in the forward scattering process, we can derive the additional differential 
photon multiplicity as follows: 

N^on.ic) = t dt p{T)ar-^\T) , (8) 

Jto 

where the cross section is given by 

= ^^f^ In , (9) 

6m dE 2m d 

with s w 2m]jE and t G (to|),s/4). Here m^ — ggT is the Debye screen 
mass. 

• Total QGP contribution to direct photons 

Combining our results, we obtain the net effect on direct 7 productions 
from the final-state jet-medium interactions in the QGP: 

D^/c{z) ^ ^ de P{e) Y^^D^/c (^3^) 
dN'' , (c) 

+ 7;-^ ' + Nl^,ic)S{l-z), (10) 
where the first term takes into account the jet quenching effect on fragmen- 

dN"' (c) 

tation photons, the second term — stands for the additional contri- 
bution from medium-induced photons which can be derived from Eq. ([7|. 
and the last term gives the correction from the inverse Compton scattering 
of the fast quark in the QGP. 

quark jet photon jet 

* (WWVW 



0000000 



Fig. 5. A fast quark jet is converted into a high energy photon jet via an interaction with a 
thermal gluon in the hot QGP. 
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4. Direct photon production in heavy-ion collisions: cold nucleeir 

matter effects 

Even though in p+A coUisions there are no final-state hot nuclear medium effects, 
these still can not be regarded as a simple superposition of p+p scatterings. There 
are different types of CNM interactions that will manifest themselves in the cross 
section for direct photon production in p+A and A+A reactions. 

• Isospin effect 

The cross sections for direct photon production for p+p, p+n and n+n 
collisions are different because they depend on the electric charges of the 
quarks (a cx e^). The different quark composition of p and n will have 
a significant impact on the nuclear modification factor of direct photon 
production R\g{pT) if the large xb region in the PDFs is probed. 

• Initial-state energy loss in cold nuclei 

A fast parton passing through cold nuclear matter may also lose energy 
before the hard scattering. There is plentiful evidence that this initial-state 
energy loss is important for heavy ion phenomenology but only recently has 
it been calculated Its effect can be modeled as: 

4>a,b/N{Xa,b, Q^) <l>a.b/N i , , ) , (H) 

where eo,fc are the fractional energy losses for the incoming partons a, b eval- 
uated in the rest frame of the corresponding target nucleus with, typically, 

ea,b < 1. 

• Cronin effect 

Initial-state multiple scattering in cold nuclear matter will broaden the 

transverse momenta of incoming partons before the hard scattering. The 

Cronin effect, in our calculation, is modelled via such kx diffusion, = 
/p \ I /t.2 \ [1] 

{k^T)med - C , (12) 

where = 0.12 GeV^, Xg — {CfICa)\ = 1 fm, L denotes the length of 
the nuclear medium and C accounts for the logarithmic dependence arising 
from large-angle scattering. Eq. ([3]) allows one to incorporate easily such 
transverse momentum broadening into the pQCD calculation. 

• Shadowing and EMC effect 

In our model shadowing was calculated from the coherent final-state par- 
ton interactions within the pQCD higher-twist coUinear factorization ap- 
proach The scale of power correction is the only parameter con- 
strained by the mean squared momentum transfer per unit length /^^/A, 
(•^2^1/3-^^^ w (2/x^L/A)g_g in minimum bias collisions, and by the world's 
data on DIS on nuclei. Only the EMC effect is parametrized in our numer- 
ical simulations. 
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Fig. 6. Nuclear modification factors K^Au^PT") ^'-'^ direct 7 production in minimum bias d+Au 
(soUd fines) and d+Cu (dashed lines) collisions at y/s = 62.4 GeV (bottom panel) and 200 GeV 
(top panel). Preliminary = 200 GeV minimum bias d+Au data is from PHENIxli^J 

5. Numerical results 

With the leading-order pQCD improved parton model in Eqs. ([1]) and ([3]), aug- 
mented by including all QGP effects, calculated in Sect. [31 and CNM effects, dis- 
cussed in Sect. 21 we can evaluate numerically the nuclear modification factor 

and compare to available experimental data at RHIC. Since we focus on hard photon 
production with > 2 — 3 GeV, the contribution from thermal photons is very 
small and thus neglected in the current study. 

Figure [6] shows numerical results of direct photon production for minimum-bias 
d-|-Au and d+Cu collisions at different colliding energies y/s = 62.4, 200 GeV by 
taking into account all relevant CNM effects, which has not yet been done be- 
fore in the literature according to our knowledge. Here to demonstrate the relative 
importance of different CNM effects, we calculated i?^^(pT) under different theoret- 
ical assumptions, namely, -R2a(?'^) with isospin and Cronin effects, Rjjiipr) with 
isospin, Cronin effect and shadowing effects, and Rlj^^ipx) with the initial-state 
energy loss as well as isospin, Cronin and shadowing effects. And in all numerical 
simulations the EMC effect is already included. We note that when px < G GeV, the 
Cronin effect is dominant and leads to cross section enhancement. For px > 6 GeV, 
the isospin effect becomes important and i?2^(pT) < 1- The effect of initial-state 
energy loss on direct 7 attenuation is substantial but the EMC effect is small and 
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noticeable only at the low colliding energy and at the largest transverse momenta. 
It is interesting to note that CNM effects have reduced direct photon production 
a,t pt ^ 15 GeV by about 25% for d+A collisions at -^/s = 200 GeV, and by about 
40% for reactions at -/s = 62.4 GeV. Due to large error bars in the data, current 
experimental measurements in d+A cannot constrain more quantitatively parton 
dynamics in large nuclei and progress in this area is urgently needed. 

Next, we show results of direct photon production in central Au+Au and Cu+Cu 
collisions at ^/s = 62.4 GeV and 200 GeV in Fig. [T] Similar to calculations for 
R'^^{pt), in Fig. [7] we demonstrate the relative strength of different nuclear matter 
effects by giving three curves of R]^^{pt)- one stands for R'J^^{pt) with coherent 
medium-induced photon bremsstrahlung and jet-photon conversion as well as initial- 
state energy loss; another curve represents _R^^(pt) with jet-photon conversion and 
initial-state energy loss; and the third curve denotes R\j^{pt) with only incoherent 
photon radiation and jet-photon conversion in the QGP. By comparing Fig. [7] to 
Fig. [6] one can see that in A+A collisions R\j^{pt) for direct photon production 
is dominated by cold nuclear matter effects, amplified by the presence of two large 
nuclei. We note that in our treatment of photon bremsstrahlung in the QGP we have 
considered the coherent interference of medium-induced photon radiation with the 
hard emission from the large scattering of the parent quark, which effect , as 
well as the initial-state energy loss effect, has been ignored in a previous study 
Experimental data clearly exclude both large Cronin enhancement and incoherent 
photon emis sion, thereby providing support for the new theoretical developments 
reported in EIH]. Copious jet-photon conversion ^ is also ruled out. In fact, in 
contrast to early speculations we find that pt < 5 GeV medium-induced photons 
contribute ~ 10% to the observed cross section and inverse Compton scattering of 
energetic quarks in the QGP is a < 25% correction. In the high pT region, the two 
enhancement contributions are very small and R\j^{pt) < 1. 

6. Conclusions 

By taking into account all relevant cold nuclear matter effects and the hot QGP 
effects we carried out the first systematic study of direct photon production in 
p+ A and A+A collisions at RHIC ^^l. After consistent modeling of initial- and final- 
state parton propagation in the medium we did not find theoretical justification for 
the previously argued large QGP-induced enhancement of the direct photon cross 
section. In fact, our results suggest that in both proton- nucleus and nucleus- nucleus 
reactions cold nuclear matter effects are dominant. Experimental data is compatible 
with such an interpretation {R\g{pT) ~ 1) and provides support for the theoretical 
underpinnings l ^ l ^ l ^^ l of our numerical results. 
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Fig. 7. R'aa^Pt) for direct photon in central Au+Au and Cu-fCu collisions at = 62.4 GeV 
(bottom panel) and y/s = 200 Ge V (t op panel) calculated including cold and hot nuclear matter 
effects. Data are taken from Refs. iMli 
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